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We analyze properties of non-hermitian matrices of size 
M constructed as square submatrices of unitary (orthogonal) 
random matrices of size N > M, distributed according to the 
Haar measure. In this way we define ensembles of random 
matrices and study the statistical properties of the spectrum 
located inside the unit circle. In the limit of large matri- 
ces, this ensemble is characterized by the ratio M /N . For 
the truncated CUE we derive analytically the joint density of 
eigenvalues and all correlation functions. In the strongly non 
unitary case universal Ginibre behaviour is found. For TV — M 
fixed and N to oo the universal resonance-width distribution 
with N — M open channels is recovered. 



I. INTRODUCTION 

Random unitary matrices may be applied to describe 
chaotic scattering |l[] , conductance in mesoscopic systems 
or statistical properties of periodically driven quan- 
tum systems (see [pi and references therein). They can 
be defined by circular ensembles of unitary matrices in- 
troduced by Dyson B. He defined circular orthogonal, 
unitary or symplectic ensembles (COE, CUE and CSE), 
which display different transformation properties pi . For 
these ensembles the distribution of matrix elements and 
their correlations are known p M . 

In the present paper we discuss properties of non- 
hermitian matrices defined as square submatrices of size 
M of unitary (orthogonal) matrices of size N, where 
N > M. These matrices may be considered as uni- 
tary (orthogonal) matrices with N — M bottom rows and 
N — M last columns truncated. Let C/rj\r,M] denote such 
a M x M matrix obtained from a unitary matrix, while 
0\n,m] is obtained by truncating an orthogonal matrix. 
The truncated matrices are non-unitary by construction, 
and their eigenvalues are located inside the unit circle. 

Motivation for such a study stems from the problems of 
chaotic scattering. Consider a mesoscopic device coupled 
to two leads, each of which supports N/2 open channels. 
The process of scattering can be described by a unitary 
S'-matrix of size N. In the diffusive regime the scattering 
matrix pertains to an appropriate circular ensemble M. 
The reflection (transmission) matrix of size M = N/2 
may be just considered as a truncation of the unitary 



S'-matrix. The random matrix approach to resonances 
in chaotic scattering was recently presented in H. In 
particular, the distribution of width of resonances in the 
presence of L open channels was derived in the weakly 
non hermitian limit for broken time reversal symmetry. 

In recent papers (H]|ll]] the authors introduce N x N 
unitary matrices enlarged in an asymmetric way to the 
size (N + L) x (N + L) by adding L upper rows and L 
last columns with all elements equal to zero. These ma- 
trices are used to describe the chaotic scattering in a ID 
model of crystal electron in ac and dc fields. It is easy to 
see, that the spectrum of such an enlarged matrix con- 
sists of 2L zeros and M — N — L complex eigenvalues 
of the truncated matrix U[n,n-l]- Our results are there- 
fore directly applicable to the problems analyzed in that 
papers. 

Related problems arise by analyzing the time evolu- 
tion of periodically perturbed systems. The model of 
the kicked rotator with absorbing boundaries was stud- 
ied in |r4-[l4[| . In this case the presence of the absorbing 
boundaries corresponds to the truncation of the infinite 
evolution matrix. Another line of research is related to 
the Frobenius-Perron operators describing the evolution 
of densities under a classical map. If the map is area pre- 
serving, the FP-operator is represented by an infinite uni- 
tary matrix, which for practical purposes is truncated to 
a finite size. Properties of the spectra of such truncated 
matrices have recently attracted a considerable attention 

mi. 

This paper is organized as follows. In section II we 
analyze truncations of orthogonal matrices and show a 
geometric interpretation of this problem. We derive the 
probability distributions of the radii of points uniformly 
covering a given hypersphere and projected into a smaller 
space. Section III is devoted to truncations of random 
unitary matrices. We demonstrate a link between the 
distributions studied and the eigenvector statistics. In 
section IV we present numerical results concerning the 
distribution of the complex eigenvalues of the truncated 
matrices . We show to what extent the ratio /i = M/N 
determines the properties of the truncated matrix. In 
section V we analytically derive the joint density of eigen- 
values for truncated matrices of CUE. From this a kernel 
is derived which determines all correlation functions. In 
particular the averaged density of eigenvalues is obtained 
for arbitrary dimensions and truncations. In the strongly 
nonunitary limit M — ► oo, N/M fixed the correlations 



of the Ginibre ensemble p3|] are obtained rescaled by the 
local mean level distance, which are thus revealed as uni- 
versal. The weakly non unitary limit N — M fixed and 
M to oo recovers the universal distribution of resonance 
widths in the weakly non-hermitian case for broken time 
reversal symmetry |9| and the corresponding correlations 
p8[ . The truncations of symmetric matrices of COE are 
briefly discussed in section VI. Convolution properties of 
the derived distributions are presented in the Appendix. 



II. SUBMATRICES OF RANDOM 
ORTHOGONAL MATRICES 

Let us start the discussion considering a simple geo- 
metric exercise. Random points cover uniformly a hy- 
persphere S^ -1 of radius 1 embedded in R N . After 
an orthogonal projection into R , where M < N, they 
are localized inside the hypersphere S M_1 or at its sur- 
face. What is the radial probability distribution Pjv,m (*)> 
where t denotes the distance of a projected point from 
the origin? 

It is helpful to analyze first the most intuitive case 
N = 3,M = 2. The surface element of the sphere S 2 
in spherical coordinates reads dfl2 = sin OdOdcj). The or- 
thogonal projection maps the points of the sphere into a 
plane. Their distance from the origin is t = sin#, which 
allows us to find the required distribution 



^3,2 (*) 



t 



2Vi~¥ 



(i) 



Analogously we get P3,i(t) = 1 for t G [0, 1] and P2,i(i) — 

l/(27rv / T^2). 

The general formula for Pn,m (t) may be obtained in a 
similar way from the element of the hypersphere S N ~ X 



df2jv-i = dipILkJi sin fc 6 k d6 k 



(2) 



Integrating out N — M variables we obtain 

r>o f,\ jJI-l/1 ,2\(N-M— 2)/2 to\ 

P N,M\ t )= c N,M t (l-i) " , (3) 

where the normalization constant can be expressed by 
the Euler beta function B(x,y) [O 



CjV,M = 



2T( 



B(f, 



N-AI 



) r(f)r( 



'N-M-- 



(4) 



Convolution relations between the distributions Pj^ M (t) 
are demonstrated in the Appendix. 

Consider an orthogonal matrix O of size N. Its first 
column can be interpreted as a vector x k = Ofci of co- 
ordinates determining a point on the hypersphere S 1 ^^ 1 . 
Let us call by Orjv,M] the upper left submatrix of O of size 
M < N. The total length of the vector represented by 

the first column of Omw and given by t = yX)fc=i x \ 
is just equal to the defined above distance of a point pro- 
jected from the hypersphere S N ~ X into the interior of 



gM-i f rom the origin. If O are distributed uniformly 
with respect to the Haar measure on O(N), than the 
points x cover uniformly the hypersphere. The distribu- 
tions P^ M (t) are then given by Eq. (J|). 

Figure l.a shows these distributions for N = 16 and 
M = 1,2,4,8 and 15. With increasing M the probabil- 
ity distribution is shifted to the right. For M = N the 
matrix remains unitary and Pi\r,iv(t) = S{t — 1). Let us 
now consider an ensemble Ow/iAn by increasing the di- 
mension N and keeping the ratio fi = M/N fixed, where 
fi < 1 . Straightforward integration allows us to compute 
for this ensemble the expectation value of t 



(t)N,/j.N 



F(f)F(^ + I) 
F(f + i)F(^)- 



(5) 



which in the limit N — > oo tends to ^JJl. The second 
moment reads (t 2 )jv,^Ar = jJi, thus the variance tends to 
zero in the limit of large matrices. This result is quite 
intuitive in view of the central limit theorem. 
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FIG. 1. Probability distribution of the radii t after projec- 
tion from R 16 to R M with M = 1,2,4,8,15. The variable 
t 2 is equal to the sum of M squared elements of a random 
orthogonal matrix distributed according to the Haar measure 
on 0(16) (a). Analogous distributions P u {t) obtained from 
random unitary matrices pertaining to (7(16), (b). 



III. SUBMATRICES OF RANDOM UNITARY 
MATRICES 

Let U[ N ^ M ] denote the MxM matrix obtained by trun- 
cation of the N x N unitary matrix U. In a similar way 

we define t — \/J2k=i l^fcil 2 - To find in this case the 
distribution P^ M (t) it is useful to represent a unitary 
matrix U(N) as a product (the normal product (x) or 
the twisted product (x)) of the hyperspheres Jl8| 



U(N) ~ S 1 x S 3 x • • • x S 2N - 3 x S 



2JV-1 



(6) 



Truncation of the dimension of a unitary matrix by one 



corresponds to the projection from S 



2JV-1 



to S 



2N-3 



which is equivalent to the truncation of the matrix 0(2N) 
by two. The same argument works for any size M of the 
truncated matrix. Therefore P$ M {t) = P% N 2 m W an d 



Pn,m(^) = C 2N,2M 



4.2 M- 



\l-t 2 ) 



N-M-l 



(7) 



with the normalization constants given by (If). Some of 
these distributions for N = 16 are plotted in Fig lb. 
Expectation values (£) are asymptotically the same for 
the both ensembles, while the variance is smaller for the 
ensemble of truncated unitary matrices E/rjv mi- 

For a fixed value of N we defined the ensembles of 
truncated matrices for each integer value of M € [1, N]. 
However, to study the evolution of a spectrum of a given 
matrix it is convenient to define an ensemble depending 
on a continuous parameter. This can be achieved in sev- 
eral different ways. For example, one may multiply the 
last column and the last vector of the matrix £/Jjv,m] by 
a parameter p g [0, 1], which mimics a continuous tran- 
sition from M to M - 1 @. 

Taking M = 1 the variable t is just the absolute value 
of the first element of a matrix \Un\. It is known Pfl9| 
that a unitary matrix of eigenvectors of a CUE matrix 
is distributed according to the Haar measure on U(N), 
while the orthogonal matrix of eigenvectors of a matrix 
typical of COE is distributed according to the Haar mea- 
sure on O(N). To establish a link with the eigenvector 
statistics let us set M = 1 and consider the distribu- 
tions Pfr A (t) and P$ t i(t) = P$ Nt2 (t). Putting y = t 2 and 
changing the variable we arrive at the known formulae 



P° N {V) 



r(f) (l-y)("-3)/2 



and 



P^(y) = (N-l)(l-y) 



JV-2 



(8) 



(9) 



which describe the eigenvector statistics for the orthogo- 
nal and the unitary ensemble |2(J,gl| . In the limit N — > oo 
they converge to the \t distributions with the number of 
degrees of freedom v equal to 1 and 2, respectively. The 
former case is often known in the literature as the Porter- 
Thomas distribution. 



IV. DISTRIBUTION OF EIGENVALUES 

Consider spectra of the truncated orthogonal matrices 
0[n,m] and truncated unitary matrices U^^m]- In both 
cases there exist M complex eigenvalues Zj = rj exp(icf>j) 
localized inside (or at) the unit circle. This is due to the 
fact that the norm of the truncation is smaller than or 
equal to the norm of the initial matrix. For the trunca- 
tions of random matrices of CUE there exist an rotational 
symmetry, U — ► C/exp(ia). Therefore P(<j>) = const, 
so we will study the radial distribution P(r). Some- 
times it is convenient to write r = e' " 2 ' and to study 
the distribution P(j) of the "level widths" 7 {?§]. For 
any fixed N the limiting cases are known: for M = 1 
the eigenvalues are trivial, r = t, so for both ensem- 
bles Pivi(r) = -Pjv,i(£)- For M = N the matrix is uni- 
tary and thus Pn Ar(r) = 5( r — 1) ol ~, m other variables, 

PnM-y) = %)• ' 





FIG. 2. Complex eigenvalues of random matrices pertain- 
ing to the ensembles U\$m (a) and E/rg 21 (b). 

Figure 2 presents 2000 eigenvales of the matrices trun- 
cated out of CUE matrices of size 5. For M — 4 there 
exist several eigenvalues close to the unit circle, while for 
stronger truncation (M = 2) the eigenvalues are clus- 
tered closer to the origin. 

In the simplest interesting case, N — 3 and M = 2 the 
data for the truncations of unitary matrices conform to 
the distribution P 3 " 2 (r) = r + 2r 3 . For comparison with 
the results of JT2J we present the numerical data as the 
distribution P{pj). The above distribution, derived in the 
following section, corresponds to the biexponential decay 



D 3" 2 (7) = ^exp(-7)+exp(-2 7 ). 



(10) 



represented by a solid line in Fig. 3a. Numerical data 
obtained for the ensemble U\$m, shown in Fig. 3b, are 
compared with the distribution ^^4(7) = ^exp(— 7) + 
\ exp(— 27) + 1 exp(— 37) +exp(— 47), which corresponds 
to P^{r) discussed below. 



P(7) 



P"(7) 




FIG. 3. Distribution of decay times P("f) for random ma- 
trices pertaining to the ensembles f7[3,2] ( a ) an d ^[5,4] (b)- 

Distributions P(r) for both ensembles obtained with 
N = 16 and some intermediate values of M are presented 



in Fig. 4. The histograms are performed out of 10 4 ran- 
dom unitary (orthogonal) matrices constructed according 
to the algorithm given in [£2| . The statistics obtained do 
not depend on which columns and rows of the initially 
unitary (orthogonal) matrix are removed during the trun- 
cation. This is due to the fact that the Haar measure 
on the unitary (or orthogonal) group is invariant with 
respect to multiplication by the permutation matrices, 
which change the order of the columns and vectors. 

With increase of M the distribution P(r) extends to 
the larger values of r. In contrast with the distributions 
P(t), for any M there exists a non-zero probability of 
finding small values of r. For small values r the dis- 
tribution P u (r) grows linearly with r. This is a purely 
geometric factor (we analyze the distribution at the com- 
plex plane), which corresponds to the uniform density of 
eigenvalues close to the origin. 
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FIG. 4. Probability distribution of the radii r of the com- 
plex eigenvalues of the M x M sub-blocks of random matrices 
of 0(16) (a). Analogous distributions P u (r) obtained from 
random unitary matrices pertaining to (7(16), (b). 

The data collected for large matrices reveal a scal- 
ing behavior: the distribution Pn,m{t) depends only on 
the ratio \i = M/N. Figure 5 shows the distributions 
P£ N / 2 i r ) an d P^ n 14(f) obtained from ensembles of ran- 



dom unitary matrices of different sizes. The larger value 
of N, the sharper is the cut-off of the probability at the 
critical radius r^ = *//Z. In analogy to the properties 
of the Ginibre ensemble one expects an infinitely sharp 
edge of the distribution in the limit N — > oo. In the case 
of large matrices the spectrum covers the entire circle of 
radius r M , while the density is largest close to the rim. 




FIG. 5. Probability distribution P u (r) of radii of eigenval- 
ues of matrices £/[jv,m] constructed from unitary matrices of 
size N = 16(o), 64(A), 256(o) and 1024(D) with M = N/4 
(a), and M — N/2 (b). Dashed lines represent the asymptotic 
cut-off at r — r M and the solid line denotes the distribution 
©■ 

For (x « 1 (and N large) the radial distribution may 
be approximated by a linear function P^ (r) ~ 2r//x with 
a cut-off at r^ . This property is characteristic of the Gini- 
bre ensemble p3j , constructed of non-Hermitian random 
matrices with no correlations between their elements. It 
is thus intuitive to expect, that for large N the con- 
straints stemming from the unitarity of U(N) do not in- 
duce very strong correlations between elements of a much 
smaller matrix of size M. 

Eigenvalues of several truncations of random orthog- 
onal matrices are shown in Fig. 6. Since the truncated 
matrix is real the eigenvalues are real or appear in com- 



plex conjugate pairs, re 4 *, re~ l( *' . Therefore these spectra 
exhibit the symmetry along the real line. Observe a clus- 
tering of eigenvalues along this line. The fraction of real 
eigenvalues equals approximatelly 0.65, 0.38 and 0.68 for 
the ensembles 0[3,2], O15M and 0[5,2]j respectively. This 
fact explains a positive probability P°[r) for r = visible 
in Fig. 4a. 





-1 1 

FIG. 6. As in Fig. 2 for truncations of random orthogonal 
matrices a) Ors 41 an d b) 0[s,2] ■ Note clustering of eigenvalues 
at the real axis. 

If the ratio fj, is kept constant, the relative number of 
real eigenvalues decreases with the matrix size. Similar 
effect is known in the theory of random polynomials. Kac 
considered the random polynomials of order M with real 
coefficients, being independent random variables drawn 
according to the normal distribution. He showed p4] 
that the fraction of real roots decreases as (\nM)/M. 
Our problem is not exactly the same since the real coeffi- 



cients of the secular polynomial of the truncated random 
matrix are not Gaussian, nor independent random vari- 
ables. In spite of this fact, our numerical results suggest 
that the fraction of the real eigenvalues of truncations 
of orthogonal matrices 0[ 2 m,ai] decreases as (ln(M)/M). 
Recent discussion of properties of random polynomials 
and their applications to quantum chaos may be found 
in p5| . The issue of clustering of zeros of random poly- 
nomials along a given curve and its relation to the time 
reversal symmetry is discussed in |26| . 

In the limit of large matrices the relative strength of 
the clustering of the complex eigenvalues along the real 
axis decreases and the distribution P(4>) becomes uni- 
form. Moreover, the radial distribution P°(r) becomes 
close to the distribution (|lj|) derived for truncations of 
unitary matrices. Although for N = 16 the differences 
between the distributions PZ(r) and P£(r) are signifi- 
cant, especially for small values of r, for large N the 
data for both ensembles seem to converge. 



V. ANALYTICAL RESULTS FOR 
TRUNCATIONS OF CUE 

Now we derive analytical results for the truncated cir- 
cular unitary ensemble. Let 



U 



A B 
C D 



be an TV x N matrix from this ensemble and A a subuni- 
tary M x M matrix. The joint density of elements of U 
can be written as 

P{U) ex S(A^A + C^C - 1)5{A^B + C f £>) x 

S(B^B + D^D-1) (11) 

with appropriate matrix 5— functions. Integrating out B 
and D we obtain as joint density of elements of A 



P(A) oc f dCS(A^A + C^C - 1) 



(12) 



with a 2M(N—M) dimensional integration over the com- 
plex parameters C. The matrix A may be brought to 
upper triangular form by a unitary transformation T : 
A = T(z + A)! 1-1 , where z is a diagonal matrix con- 
sisting of the complex eigenvalues of A and A is strictly 
upper triangular ( Schur decomposition). The transfor- 
mation can be made unique restricting T to a certain 
cosetspace. The Jacobian of this transformation is g iven 
by the square of the Vandermonde determinant ||27| 



W\ 






such that after integrating out the unitary transforma- 
tions T the joint density of eigenvalues is given by 



P{z) oc \V\ 2 dC dA<5((z t + A f )(z + A) + &C - 1). 

(13) 

In the following we first integrate out the M(M — l)/2 
complex parameters Ay and then the M(N — M) com- 
plex parameters C written as complex vectors CV For 
i < j we have the hierarchical equations 



^A tJ +C}C J+ J2^^=0. 



(14) 



k<i 



Integration over A yields the Jacobian Ili<j \ z i\ 2 an d a 



product of <5— functions 

HSdzif + CJXid-l). 



(15) 



where JQ denotes an N — M x N — M matrix defined by 
the quadratic form CJXiCi which is given by 

k<i 

containing the solution A^^ of equ.(^4J) and depending 
otherwise on C'k only for k < i. The integration over 
d can now be done successively starting from Cm and 
yields the factor 



M 



IB 1 "!*! 2 ) 



2\N-M-1 



9(l-|z J | 2 )/det(X l ) 



where 6(.) denotes the Heaviside step function, <d(x) = 1 
for x > and zero otherwise. 

For det(Xj) we can derive from the equations ( O ) for 
A and using implicitly the S— functions (Eq) the recursive 
relation 

dct{X l+1 ) = dct(X i )/\z l \ 2 

with det(Xi) = 1. Thus the previous Jacobian 
Til |.Zi|~ 2 will be compensated and the final very simple 
and important result is 

1..M M 



p{z)k n^-^nc 1 -!*! 2 ) 



2\N-M-\ 



i<3 



e(i-N 2 ) 



(16) 



This result is completely analogous to the Ginibre ensem- 
ble and we know immediately all correlation functions by 
the method of orthogonal polynomials g . Here the pow- 
ers z" _1 arc already orthogonal. An equivalent method is 
to consider the joint density P(z) as the absolute square 
of a Slater determinant of normalized wave functions 



with 



<j )n {z)^z n - 1 w{\z\ 2 )/^N~ 

,(x) = (i-x)( N - M - 1 v 2 e(i-x), 



where N n stands for an normalization factor. 

The kernel, which determines all correlation functions 
is i0 

M 

K( Zl ,z* 2 ) = ^(z 1 z*)^ 1 W (|z 1 | 2 ) W (|z 2 | 2 )/iV„. 

n=l 

For example the cluster function is given by Y{z\, z<i) — 
\K(zi,Z2)\ 2 and the averaged density of eigenvalues z 
normalized to 1 is given by 

M 

p(z) = K(z 7 z*)/M = - £ \z\ 2n - 2 w 2 (\z\ 2 )/N n . (17) 

n=l 

The normalisation factor N n is easily calculated as 

N n = n(n - 1)!(JV -M- 1)\/{N -M + n-l)\ 

For example with r 2 = \z\ 2 we obtain for the distribution 
of r with M = N - 1 



P(r) = — (r + 2r 3 - 
v ; AT 

and in general with x = r 2 
2r (1-x) 



3r b 



Mr 



2M-1 



^(r) 



iV-M-l / J \ W-M n _ ^jy^ 



M (iV-M-l)! VA 



(1 - X N ) 

l-x 



(18) 



There are two important limiting cases for large M: 
H = M/N fixed and L = N - M fixed. For fixed p and 
M to oo we find the mentioned scaling behaviour: 



P(r) = (--!) — 

K ' > j (l-r 2 ) 



2 "1 2 



(19) 



for r 2 < p and P(r) = otherwise. The distribution 
shows a gap near the unit circle. This gap resembles 
the one obtained for resonances in the chaotic scatter- 
ing problem for large number of channels (£9). In this 
strongly non unitary limit we are also able to simplify 
the cluster function 

Y(z, z + S) = {Mp{z)f exp(-TrMp(z) \6\ 2 ) (20) 

which is just the Ginibre behaviour P3lp| with the 
dist ance 6 rescaled by the local mean level distance 
l/y/Mp(z) given by equ.(p^) through p(z) — P(r)/2irr. 
The same can be shown for the nearest neighbour dis- 
tance distribution obtained by Grobe et al. |31[ and ap- 
plied to a damped chaotic kicked top. 

In the other limit of fixed L = N — M and M to oo, 
which may be considered as weakly nonunitary, we re- 
cover exactly the universal resonance-width distribution 
M for perfect coupling to L channels with y = N(l — r) 



Similarly the cluster function obtained in this limit can 
be shown to coincide with the one obtained by Fyodorov 
and Khorushenko Eq] for chaotic scattering with a finite 
number of perfectly coupled channels. The statistics (|2l| ) 
has also been found by Kottos and Smilansky [p2| for 
chaotic scattering on graphs and by Gluck et al. RIO 11 
for a model of crystal electron in the presence of dc and 
ac fields. In both of these works the S-matrix is reduced 
to the resolvent of a subunitary matrix as is investigated 
in the present paper. 



VI. SUBMATRICES OF UNITARY SYMMETRIC 
MATRICES 

For several applications one uses symmetric unitary 
matrices typical for the circular orthogonal ensemble. 
This case is relevant if the physical system possesses time 
reversal symmetry, or any generalized anti-unitary sym- 
metry |3|]. Let U be a random unitary matrix typical 
of CUE. It is easy to prove that the symmetric matrix 
W := UU T is typical to COE @j. We shall thus define 
the ensemble of truncated symmetric unitary matrices 
W[jv,Af] ■ I* 1 the definition of this ensemble the position 
of the submatrix is crucial. We take the left upper part 
of the symmetric matrix W , thus the truncated matrices 
WwM] are symmetric. 

The distributions P(t) and P(r) for the symmetric ma- 
trices generated out of COE matrices of size N = 16 are 
shown in Fig. 7. Each plot contains data from 10 sym- 
metric random unitary matrices. Note the differences 
between these figures and the corresponding data for or- 
thogonal and unitary matrices presented in figures 1 and 
4. If the truncation of the matrix W is performed asym- 
metrically (e.g. we take the lower left submatrix), the 
distribution P(r) becomes closer to this corresponding 
to the truncations of random unitary matrices E/rjv,Ml- 

In the asymptotic limit the properties of the ensemble 
of the truncations of symmetric matrices W[n,m] depends 
on the same scaling parameter p = M/N. Moreover, the 
distribution P(r) becomes close to the corresponding one 
for the unitary ensemble described by the distribution 
(J19|). Therefore we may conjecture that the distribu- 
tion P M (r) , which describes the distribution of moduli of 
eigenvalues of truncated matrices in the limit of large TV, 
is universal and does not depend on the initial ensemble 
of random matrices, provided M/N is fixed. This corre- 
sponds to results for resonances in the limit of L/N fixed 
and TV to oo |g9[ . In the contrary, there are differences to 
be expected in the limit of weakly nonunitary matrices: 
N-M fixed and M to oo fJoLEl. 
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FIG. 7. Probability distributions P(i) (a) and P(r) (b) for 
complex symmetric matrices W[i6,a/] with M — 2,4, 8, 15. 



VII. CONCLUDING REMARKS 

Three families of ensembles of random matrices are 
proposed. They are defined by cutting an M dimensional 
submatrix of an initially ./V-dimensional unitary matrix, 
pertaining to a given ensemble of unitary, unitary sym- 
metric or orthogonal matrices. Using a link between the 
truncation of an orthogonal matrix and the projection of 
a hypersphere into a smaller dimensional space we found 
the probability distributions of the lengths t of the pro- 
jected vectors. 

Truncated matrices are not unitary and their complex 
eigenvalues are located inside the unit circle. We derived 
an analytical formula for the distribution P(r) of mod- 
uli of eigenvalues of truncations of the CUE matrices. It 
takes a particularly simple form for small values of N 
and M. In the asymptotic limit N — > oo this distribu- 
tion depends only on the scaling parameter fi — M/N, 
provided M/N is not very close to 1. For small r the 



distribution P^ (r) grows linearly, later displays a nonlin- 
ear behaviour and eventually suffers a sudden cut-off at 
r fj, = \JV- F° r N » 1 the probability distribution P(r) 
does not depend, whether the initial matrices are orthog- 
onal, unitary or unitary symmetric, again if M/N is not 
very close to 1. In this strongly nonunitary limit at least 
for the case of broken time reversal symmetry correla- 
tions are shown to coincide with those obtained from the 
Ginibre ensemble of general complex matrices, provided 
distances are rescaled by the local mean level distance. 
In the weakly nonunitary limit N — M fixed and M to oo 
the eigenvalue distribution again for broken time reversal 
symmetry is shown to coincide with the universal reso- 
nance widths distribution in the weakly non-hermitian 
limit. The same is true for the 2-point cluster function. 
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APPENDIX A: CONVOLUTION PROPERTIES 
OF THE DISTRIBUTIONS P{t) 

In the appendix we demonstrate the convolution prop- 
erties of the distributions Pjy M (t) which might be used 
to derive the formula (||). We start thus with a random 
orthogonal matrix O(N) or with random points covering 
uniformly the hypersphere S' JV_1 of radius 1. Their dis- 
tribution in the polar coordinates is given by Eq. (|2j). 
For simplicity we will denote the distance from the ori- 
gin of a point projected into R M by ijv,M- It is just 
the argument of the distribution P^ M (t). Due to the 
definition of the polar coordinates ijv,JV— l — sin#7v-2, 
tN,N-2 = sin0jv-2sin0jv-3, • • •, and ijv,i = cos9 N ^ 2 - 
Therefore all variables i^A/ may be rewritten as the 
product consisting of L = N — M factors 



N-M 



tN,M — 1 [ ^JV- 



k+l,N-k- 



(Al) 



fc=i 



This factorization allows us to find the distributions (0). 
Probability distribution of a sum of two independent 
random variables z = x + y is given by the standard 
convolution P x o P y := P{z) = J_ oc P x (x)P y (z — x)dx. 
In a similar way, the distribution of the product of two 
independent random variables, z = xy, is given by the 
product convolution 



P x * P y := P(z) = 



P x (x)P y (-)^-dx. (A2) 

v ar \x\ 



In the general case the integration should be performed 
over the entire real axis, but in our case the intergration 
is restricted to \z , 11 . since all arguments t £ [0,1]. 

Factorization (Al) allows us to write convolution re- 
lations between probability distributions P^m^)- For 
example 



°3lW — °32 *P; 21 



"42 W ~ P 43 * "32 ) 



(A3) 



(A4) 



*?i (*) = ■*$» **&**&• (A5) 

In general we obtain a convolution relation 

Pn,mW — Pn,n~i * Pn-i,n-2 * ' ' ' *°m+i,m 1 (-"-6) 
which might be used to derive formula (H). 
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